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Abstract. The effective potential between charged colloids trapped at water 
interfaces is analyzed. It consists of a repulsive electrostatic and an attractive capillary 
part which asymptotically both show dipole-like behavior. For sufficiently large colloid 
charges, the capillary attraction dominates at large separations. The total effective 
potential exhibits a minimum at intermediate separations if the Debye screening length 
of water and the colloid radius are of comparable size. 
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In view of various basic and applied issues such as the study of two-dimensional melting 
PP, investigations of mesoscale structure formation [2 j or engineering of colloidal crystals 
on spherical surfaces [3], the self-assembly of sub-/zm colloidal particles at water-air 
or water-oil interfaces has gained much interest in recent years. These particles are 
trapped at the interface if water wets the colloid only partially. This configuration is 
stable against thermal fluctuations. It appears to be even the global equilibrium state, 
because it is observed experimentally that the colloids immersed in the bulk phases 
are attracted towards the interface 1J. For charge-stabilized colloids at interfaces, the 
repulsive part of their mutual interaction is well understood and resembles a dipole- 
dipole interaction at large separations. This asymptotic interaction is caused by charges 
at the colloid-water interface [3] or by isolated charges at the colloid-air (or oil) interface 
Nonetheless, charged colloids at interfaces exhibit also attractions far beyond the 
range of van-der-Waals forces. According to Refs. (fU HI El El EH IH] , polystyrene spheres 
(radii R = 0.25 . . . 2.5 /im) on flat water-air interfaces using deionized water exhibit 
spontaneous formation of complicated metastable mesostructures. They are consistent 
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with the presence of an attractive, secondary minimum in the effective intercolloidal 
potential at separations d/R ~ 3 ... 20 with a depth of a few k%T . In Ref. ^2]; PMMA 
spheres with radius R = 0.75 /im were trapped at the interface of water droplets 
immersed in oil. Here, the secondary minimum has been detected at a separation 
d/R = 7.6 and is reported to be surprisingly steep. The tentative explanation of these 
findings given in Ref. invokes an analogue of long-ranged flotation or capilllary 
forces which decay oc 1/d. This interpretation was criticized in Refs. [TB*] HI] which 
both concluded that possible capillary forces in this system are much shorter ranged, 
i.e., oc d~ 7 , but the authors of these references disagree with respect to the sign of 
these shorter-ranged forces. Recently we have shown ^Hj quite generally that long- 
ranged flotation-like forces oc 1/d can only arise in mechanically non-isolated systems. 
For isolated systems the capillary force is indeed much shorter-ranged and within a 
superposition approximation the power law discussed in Ref. |Hj is found. For the 
experiment involving the interface of a mesoscopic droplet jT^j mechanical isolation 
may indeed be weakly violated and thus small flotation-like forces can appear. Their 
interplay with capillary forces arising from the droplet curvature is not clear yet and 
is currently under scrutiny jTH]- On the other hand, for the experiments performed 
on flat interfaces [HJ 13 El M EDI HI] , mechanical isolation holds if there is no external 
electric field present and thus for them flotation-like forces cannot give rise to the 
observed attractions. Out of other attempts to explain the nature of this colloidal 
pattern formation we mention Ref. J7j, in which this is attributed to oil contaminations 
of the interface, and Ref. [TU] . in which colloidal roughness is proposed as a source of 
attractive capillary interactions. At present, these attempts to explain the observed 
colloidal patterns at interfaces are only of qualitative nature. 

Thus a theoretically sound mechanism for the appearance of an attractive minimum 
in the intercolloidal potential at large separations has not been found yet. Here we 
analyze the interaction between colloids at interfaces within the approach developed 
in Ref. for a mechanically isolated system and we provide conditions for an 
asymptotically attractive intercolloidal potential and for the appearance of such a 
minimum. 

In going beyond the superposition approximation studied in Ref. [IK] , we derive 
the full capillary interaction potential between two colloids as a functional of a general 
stress field acting on the interface. This capillary potential is studied for two cases: (i) 
the Debye screening length of water, is much smaller than the colloid radius R, 
and (ii) R and k~ 1 are of comparable size. For a sufficiently high charge on the colloids, 
in both cases the ensuing capillary attraction turns out to be asymptotically stronger 
than the direct repulsion. Moreover, in case (ii) a minimum in the total (repulsive plus 
capillary) intercolloidal potential is found at intermediate separation Krf min > 10. 

We consider two spherical colloids a = 1,2 trapped at a deformed interface 
(meniscus) with vertical coordinate z = u(r = (x,y)). We denote h a as the vertical 
coordinate of the center of colloid a and r a its lateral position so that d = \r± — T2\. 
We define a reference configuration (with respect to which free energy differences are 
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measured) by a flat interface u = and h a re{ = —Rcos9. Here, 9 is Young's angle 
and thus in the reference configuration the colloids are vertically positioned such that 
Young's law holds at the horizontal three-phase contact circle with radius ro, r ef = -Rsin 9. 
The corresponding free energy is given by 

f = 1 \Sr ^l + ^-^u + Yilir—t d£[Ah a -u} 2 -F a Ah a \. (1) 



The first line of Eq. (0) comprises free energy differences associated with the change 
in meniscus area, in meniscus gravitational energy (7 is the water-air surface tension 
and A is the capillary length) and with forces on the meniscus; the stress IT denotes 
the vertical force per unit area on the meniscus in the reference configuration. The 
first term in the second line takes into account the changes in water-colloid and air- 
colloid surface energies; Ah a is the difference in colloid center position with respect 
to the reference configuration. The second term describes the energy difference if 
colloid a is shifted by the force F a (which is evaluated in the reference configuration). 
S a is the circular disk of radius r 0ire f delimited by the three-phase contact line dS a 
formed on the colloid a by the fluid interface in the reference configuration (traced 
counterclockwise), and S men = M 2 \(Si \J S 2 ). This expression for the free energy is valid 
as long as the deviations from the reference configuration are small: |w|/Vo ire f, | Vw| <C 1. 
A sufficient condition for this is \ep\ <C 1, with ep := —F a j 27r7ro, re f. Mechanical 
isolation implies 2F a = — [ Q d 2 rll, i.e., the forces on the colloids are balanced by 
the force on the meniscus [15J. The meniscus-induced effective potential between the 
colloids is defined as V men (d) = F(u, h a \ eq ] d) — T{u, h a \ eq ; d — > 00). The equilibrium 
free energy is found by minimizing Eq. ((TJ with respect to u and h a . The equilibrium 
meniscus shape for the two colloids being infinitely apart is given by the superposition 
u(d — > 00) = m + U2 of the single colloid menisci, obtained in the presence of the stress 
field fl(d — > 00) = ni+n 2 . Here, u a = u(\r— r a \) is the equilibrium meniscus around one 
colloid and likewise U a = n(|r — r Q |) is the stress field on the interface caused by a single 
colloid. Furthermore we define the single-colloid quantity Bp = f M2 \ S d 2 rIL/ 27r7r 0ire f. 
For finite d, the stress field can be written as II(g?) = Hi + U 2 + 2 n m and likewise we 
decompose u(d) = U\ + u 2 + u m . Minimization of Eq. (JTJ) yields 

V u m - — = , (2) 

A z 7 

with the boundary condition u m = for r — > 00, and at r G dS a {(3 ^ a): 
d(u m + U(3) u m + up - (u m + Up) a 

o — £p-£F-\ , l-Jj 

on a r ,ref 

with 27rr ref (-) a = <f ds (-)d£. The superposition approximation entails n m = u m = 
for all d (i.e., ep = ep), and was analyzed in Ref. To identify the corrections to 

it, we write Vm en = V sup [n] + V^pTj and obtain after some algebra (in the following we 
consider A — ► 00; this limit can be safely taken in the case of mechanical isolation [To] ) 



V sup = / d 2 r n 2 M 2 - d 2 r Tl^ + 27T7r , rc f e F (u 2 )i 
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V m = - jd 2 r (Tl 2 u m + 2 Tl m u 2 + U m u m ) + rei (4) 

" 5'mcn 



x 



,2 J Uf Ml , , (U 1 + U 2 ) 1 

[e F - Ep) + ep [e F - e F ] 



T O.ref ' T o,ref 

In the following, we apply this general expression for the capillary potential within 
two electrostatic models which provide explicit expressions for IT and Il m . Generically, 
the water phase contains screening ions which lead to a finite Debye screening length 
n~ l = [e 2 j 87r/3coe 2 ) 1//2 where Cq is the concentration of monovalent ions, e is the 
elementary charge, and = ksT. We denote by ex, e 2 , e c the permittivities of air, 
water, and the colloid, respectively, using Gaufi units. 

1. nr 0)Te i 1. The electrostatics of the single-colloid configuration has been 
analyzed in Ref. [H] with the result that only charges at the colloid-air interface generate 
a stress on the interface. The potential <T>o a t the air-water interface is small so that 
the interface resembles a perfect conductor. According to Ref. [TSj . IT ~ (ex/ 8n)E 2 0+ , 
where E Z:0 + denotes the ^-component of the electric field right above the interface, 
leading to the following robust parametrization of the corresponding numerical results: 

m = 3^)(j--iY 1 (j-)' 1 ^ , (5) 

?~0,ref V r 0,rcf / \ r 0,ref/ 

where ep > 0, b = + l)(fi + 2)(/i + 3)/6 and \i G (0,1) is a fitting parameter 
depending on e c and the contact angle 6. IT has an integrable divergence as r — > r re{ 
and rapidly (i.e., r > 2r 0ire f) reaches its asymptotic dipole behavior cx r -6 . The dipole 
asymptotics renders the repulsive part of the intercolloidal potential V^ ep for large d, 

Vre P = ^rl e ,e F b^)^f) 3 , (6) 

and IT m = (n 1 n 2 ) 1 ' /2 because the electric field of the two-colloid configuration is the 
superposition of the electric fields in the single-colloid configurations. This stress is 
strongly peaked around the colloid centers and the main contribution in Eqs. (JH) stems 
from the regions around the colloids as d — > oo. Therefore, in order to obtain V men to 
leading order in 1/d one can employ the approximation IT m w IT 1 / 2 ((i)(IT^ 2 + IT^ 2 ) to 
solve Eqs. (0) and Q and to evaluate the integrals in Eq. (jlj): 

V men (d) « -2n ir l e{ el (^f) 3 ^ [1 + M{n)\ . (7) 

Here, the function M(fi) is given by a lengthy analytical expression; it increases almost 
linearly for \x G (0,1) with M(0) = and M(l) = 1/5. The leading asymptotic 
behavior oc d~ 3 stems from V m (the superposition approximation predicts l^ up oc d~ 6 
|15p. Repulsive and capillary force decay with the same power but with opposite 
signs of the amplitudes. Hence, the total intercolloidal potential will be attractive if 
£ F > e FiCrit = (ji + 1)/4[1 + M(/i)] with 1/4 < s F}Clit < 5/12 for < // < 1. The 
appearance of capillary attractions thus depends sensitively on the colloidal charge 
via ep. These critical values for ep are at the limit of applicability of this model: 
The maximum gradient and the maximum deviation for the meniscus occur at the 
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three-phase contact line and are given to leading order in 1/d by the single-colloid 
solution: |Vw| ~ £p, u/r ^ TC f m —fiEp/4. Experiments which are performed in the limit 
o,ref > 1 H CHj estimate charge densities on the colloid-air surface compatible with 
< Ep <C 1 and thus the electrostatic repulsion would always be stronger than the 
capillary attraction. Even if ef were large enough such that the capillary attraction 
would dominate, Eqs. © and (JJJ) would not render a minimum in the total potential at 
intermediate separations. However, the possibility for the occurence of such a minimum 
arises outside the regime M" 0ire f 3> 1- 

2. M" 0)ro f < 1. In this regime, n~ l provides an additional length scale which leads 
to interesting crossover phenomena and the charges on the colloid-water surface are not 
completely screened in the range of separations of interest. The stress is given by 



Here, E z $+ and Euq are the perpendicular and the parallel component of the electric 
field at the interface on the air side. The osmotic pressure p osm = j3~ l An is generated 
by the excess ion number density An at the interface on the water side. In order to 
solve the electrostatic problem, we introduce two simplifications: the Debye-Hiickel 
approximation, and the point-charge approximation, by which the total charge q 
of the colloid is concentrated at the center of S a . Then the osmotic pressure is 
given by p osm = (£2/ 87t)k 2 $q and the electrostatic repulsion between two colloids by 
Vre P (d) = q&o (d), where <3> is the electrostatic potential at the interface. V Tep exhibits 
a crossover from a screened Coulomb repulsion to a dipole repulsion at nd c ~ 8 ... 10 
The single-colloid stress, Eq. (JEJ), is dominated by Eu and p osm for nr < 6 . . . 8, 
yielding Tl(r) oc exp(— 2/tr)/r 4 . For larger distances, it is dominated by E z q+ and the 
familiar dipole form arises, U(r) oc r~ 6 . The stress II exerted by two colloids can again 
be determined by superimposing the electric fields and potentials of the single-colloid 
solution, leading to 



n - = i 1 1 - - (£*,o+)i(£*,o+)2 + -V -(^11,0)1(^11,0)2 COS0+ ^ 2 ($o)i($o) 2 -(9) 



As before, the subscript a = 1, 2 denotes evaluation of the single-colloid function at 
|r — r a | and cos0 = (r — r x ) • (r — r 2 )/|r — r x | |r — r 2 |. In the limit d 3> d c , u m and V^en 
can be estimated via an expansion in 1/d using a peak approximation for n m as before. 
This yields an asymptotically attractive capillary potential which again decays oc d~ 3 . 
For intermediate separations d, we have compared these results with a full numerical 
solution to Eqs. (J2J) and Q without any approximation for n m or V men . In Fig. d (left 
panel) we show the results for Vmen: For d > d c the l/<i-expansion and the numerical 
results agree with each other, both exhibiting the dipole-like asymptotic behavior. For 
smaller d the extrapolated 1 /(^-expansion predicts a repulsive capillary potential which 
is borne out by the numerical results only partly as Kr 0ire f — > 1. The superposition 
contribution V snp is again negligible in the expression for Vmen- The total intercolloidal 
potential V tot = V men + V rep is depicted in Fig. [T] (right panel) for the exemplary value 




(8) 
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Figure 1. (colour online) The capillary potential (left panel) obtained from the full 
numerical solution in units of Vo = q 2 K 2 r q jle { / 2tt 62 for Kra iTe { = 0.1 . . . 1.0 and the 
parameters £2/61 = 81 and k\ = 100. The peak approximation fails for nd < 8. The 
total intcrcolloidal potential Vtot — V mcn + V lep (right panel) for ep = 0.6 and the same 
parameters and units. 




e F = 0.6. Asymptotically V tot is attractive for e F > e FjC rs 0.34[1 + 0.30/(/tr 0>rof ) 2 ] 1 / 2 , 
obtained from a fit to numerical results for nr 0>Te f £ [0.1, 2.0]. 

A quantitative comparison with experiments is difficult because the three important 
quantities q, 8, and k , which enter into e F and the scale of Vmen have not been 
determined separately for the same system. Therefore, we can only estimate whether 
a minimum as obtained in Fig. ^ can occur in actual experiments. For almost all 
experiments pure water is claimed to have been used, for which n^ 1 « 1 /mi. The 
total charge is given by q = 2ira R 2 (1 + cos 6) where a is the charge density. Thus, 
for colloids with R = 0.5 /mi and 9 = it /2 at pure water-air interfaces (7 = 0.07 
N/m), we obtain numerically e F m 45(<jnm 2 /e) 2 . The potential scale at T = 300 
K is given by V Q = g 2 K 2 r , re f/(27re 2 ) « 1.4 ■ 10 8 k B T {a nm 2 /e) 2 . The value e F = 0.6 
used in Fig. ^ corresponds to a charge density a = 0.12 e/nm 2 (literature estimates 
for the actual charge density vary between 0.07 jH] and 0.53 ^Uj e/nm 2 ), in which case 
Vo = 2-10 6 ksT. For M" ,ref = 0.5 the minimum in the total potential (see Fig.QJ) occurs 
at nd min m 13 with V m i n ~ —1.6 • 10 Vq = 32 k B T. This minimum is shallow enough 
that thermal movements of the colloids around the minimum position should be visible, 
similar to reports in the literature. 

Such large charge densities actually call for solving the Poisson-Boltzmann equation 
near the colloids. Incorporating the correct geometry, this is a very involved numerical 
task and will be considered in future work. From that we expect the following 
modifications: The screening ions will concentrate near the colloid in a layer with 
thickness of a few nm (as given by the Gouy-Chapman length Iq = 81162/ (flea)). Outside 
this layer, the colloids appear as heavily screened objects with an effective charge q e s and 
the Debye-Huckel approximation is valid. This charge g e fr enters both the electrostatic 
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repulsion and the capillary attraction such that asymptotically V Tep oc ql s (r 0tTe {/d) 3 and 
V men oc Sp g 2 ff (r 0jre f/<f) 3 . The total force on one colloid (as expressed by ep) is not 
determined by q e g, but rather by the interfacial stress in the screening layer (dominated 
by p os and -E^o) • Since this layer is thin compared to the colloid radius, near the three- 
phase contact line the problem appears to be similar to that of a charged plate half- 
immersed in water. For an order of magnitude estimate of ep we have used the analytical 
solution for a fully immersed plate and obtain ep ~ e 2 / ' {^^l3 2 e 2 lo) = 1.2(<rnm 2 /e) for 
the parameters discussed above. Thus, Ef = 0.6 (as used in Fig. [T| is obtained for 
a = 0.5e/nm 2 , not too far from the Debye-Hiickel result. We emphasize that the 
occurrence of a potential minimum for ep > ep >c (i.e., for a sufficiently large charge 
q) is a consequence of the intermediate distance crossover in II from being pressure- 
dominated (by i^i^o and p OS m) to tension-dominated (by E z ^+). This mechanism is 
captured correctly by the presently employed approximations (point charge and Debye- 
Hiickel treatment). 

To summarize, we have calculated the effective intercolloidal potential for charged 
colloids floating on a water interface. We have derived a general expression for the 
capillary potential as an explicit functional of the stress on the interface. We have 
quantitatively studied the capillary potential when the stress is due to the electrostatic 
field of charged colloids in the cases that the radius R of the colloid compared with the 
Debye screening length of water is either very large or of about the same size. In both 
cases the asymptotic behavior of the capillary potential and of the direct electrostatic 
repulsion is equal, oc d~ 3 , but different in sign. The superposition approximation, 
predicting a capillary potential oc d~ 6 , is insufficient because it takes into account 
only the energy change as the subsytem "one colloid+surrounding meniscus" is shifted 
vertically in its own (single-colloid) force field. The correction to the superposition 
approximation, embodied in II m , considers the additional force by the second colloid. 
The capillary attraction only dominates for a sufficiently large charge and the total 
intercolloidal potential exhibits an attractive minimum only if k~ 1 ~ R. This minimum 
can be understood as a consequence of the pressure-to-tension crossover in the stress 
acting on the interface. The depth of the potential minimum is significantly reduced 
compared to the natural energy scale 7.R 2 of capillary interactions and is of the order 
of several &bT for parameters relevant for actual experimental conditions. 
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